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ABSTRACT 

We investigate the physics of gas accretion in young steUar clusters. Accretion in 
clusters is a dynamic phenomenon as both the stars and the gas respond to the same 
gravitational potential. Accretion rates are highly non-uniform with stars nearer the 
centre of the cluster, where gas densities are higher, accreting more than others. This 
competitive accretion naturally results in both initial mass segregation and a spectrum 
of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal- 
lobe radius instead of the commonly used Bondi-Hoyle accretion radius. This works 
as both the stellar and gas velocities are under the influence of the same gravitational 
potential and are thus comparable. The low relative velocity that results means that 
-Rtidai < Rbh in these systems. In contrast, when the stars dominate the potential and 
are virialised, Rbh < -Rtidai and Bondi-Hoyle accretion is a better fit to the accretion 
rates. 
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1 INTRODUCTION 

Star formation involves the collapse of molecular cloud cores 
under their self-gravity through >20 orders of magnitude in 
density. This process is extremely non-homologous in that 
a small fraction reaches stellar densities while the vast ma- 
jority of the mass is still infalling (Larson 1969). Thus, a 
stellar core, comprising a small fraction of a solar mass, is 
initially formed which grows through the accretion of the 
infalling envelope. It is the accretion which ultimately deter- 
mines the final stellar properties (Stahler, Shu & Taam 1980; 
Palla 1999) as well as the observed characteristics during 
its pre-main sequence evolution (Henricksen Andre & Bon- 
temps 1997; Andre, Ward- Thompson & Barsony 2000) . 

Complicating this picture is the fact that most stars do 
not form in isolation but rather in groups from binary sys- 
tems (Mathieu et al. 2000) to stellar clusters (Clarke, Bon- 
nell & Hillenbrand 2000). The separation between the indi- 
vidual components of these systems is typically less than the 
size of an accreting envelope such that they must compete 
for the reservoir of material. 

Surveys of star forming regions have found that the ma- 
jority of pre-main sequence stars (50 to 90 per cent depend- 
ing on the region considered) are found in clusters (e.g Lada 
et. al. 1991; Lada, Strom & Myers 1993). These clusters con- 
tain anywhere from tens to thousands of stars with typical 
numbers of around a hundred (Lada et. al. 1991; Phelps & 
Lada 1997; Clarke et. al. 2000). Studies of the stellar content 
of young (ages « 10*^ years) clusters (e.g. Hillenbrand 1997) 



reveal that they contain both low and high-mass stars in 
a similar proportion to that in a field-star IMF (Hillen- 
brand 1997). Observations of young clusters have shown that 
they are usually associated with massive clumps of molec- 
ular gas (Lada 1992). Typically, the mass of the gas in the 
youngest clusters is greater than that in stars (Lada 1991), 
with up to 90 % of the cluster mass in the form of gas. Gas 
can thus play a crucial role in the evolution of the clusters. 
Furthermore, there is a degree of mass segregation present 
in the clusters with the most massive stars generally found 
in the cluster cores (Hillenbrand & Hartmann 1998; Car- 
penter et al. 1997). These systems are generally too young 
for two-body relaxation to explain these observations im- 
plying that mass segregation is an initial condition of stellar 
clusters (Bonnell & Davies 1998). 

Previous studies have explored the physics of accretion 
in binary systems (Bate & Bonnell 1997; Bate 2000) and 
in small stellar clusters (Bonnell et al. 1997). Accretion in 
binary systems plays an important role in determining the 
system's mass ratio and separation (Bate 2000) . In the more 
chaotic environment of a small stellar cluster, the accretion 
rates are highly non-uniform with a few stars accreting much 
more than the rest. The accretion rates primarily depend on 
the star's position within the cluster such that those near- 
est the centre have the highest accretion rates (Bonnell et 
al. 1997). This competitive accretion results in a spectrum 
of stellar masses, from initially equal masses, and as such 
is a promising candidate for explaining the observed initial 
mass function (Zinnecker 1982; Bonnell 2000). 
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In this paper we are extending the work on smaUer stel- 
lar clusters to larger systems containing up to one hundred 
stars. The aim is to explore the physics of the competitive 
accretion in order to get a better grasp of the likely effects 
of competitive accretion on the spectrum of stellar masses. 
In § 2 we discuss the numerical calculations. In § 3 we dis- 
cuss how the accretion is related to the cluster dynamics. 
Section 4 presents the results on accretion and mass segre- 
gation while § 5 presents our results on modelling the physics 
of competitive accretion. Our conclusions are presented in 
§ 6. 



2 CALCULATIONS 

The calculations presented here were performed with a hy- 
brid SPH-Nbody code to model the combined presence 
of stars and gas and their mutual interactions. This code 
is based on a standard 3-D Smoothed Particle Hydrody- 
namics (SPH) code (Benz 1990) with gravitational inter- 
actions calculated via a tree-code (Benz et al. 1990). The 
gas is modelled by the normal SPH particles while the 
stars are modelled using sink-particles (Bate, Bonnell & 
Price 1995). These sink-particles interact with other par- 
ticles only though gravitational forces and the accretion of 
the gaseous SPH particles. The accretion is modelled by re- 
moving gas particles within a predetermined radius of one 
of the sink-particles, providing they are bound to the sink- 
particle, and adding their mass and momentum to the sink- 
particle's (Bate et al. 1995). The radius at which this accre- 
tion occurs, the sink-radius -Rsink is chosen to be sufficiently 
small that it does not overtly affect the gas flow outside 
this radius. In practice this means that it must be smaller 
than any physical accretion radius (e.g. Bondi-Hoyle accre- 
tion radius) where the gas becomes bound to the star. The 
sink-radius cannot be too small either as this translates into 
a prohibitively small time-step. In this study, we use a sink- 
radius of -Rsink ~ 10~^ X -Rciust, the initial cluster radius. 

The simulations were performed with 10^ SPH parti- 
cles to model the gas. The end-point of the simulations was 
generally chosen to be when roughly half of the SPH parti- 
cles were accreted, although several were run till all the gas 
was accreted. The analysis was performed before the ma- 
jority of the gas was accreted in order to ensure sufficient 
resolution. Furthermore, several simulations were rerun with 
accreted particles each replaced by nine lower-mass particles 
to ensure that the results were not overtly affected by the 
resolution. This was indeed found to be the case in all the 
simulations except those with both 'hot' gas and stars where 
the low accretion rates require larger particle numbers to be 
fully resolved. Limitations of the simulations presented here 
are that there is no feedback from the stars onto the gas 
and that the gas does not initially contain any turbulent (or 
other) motions. These two effects could play a significant 
role in the cluster evolution and in the accretion process but 
are beyond the scope of this paper. 



2.1 Initial Conditions 

We have simulated accretion onto a large number of clus- 
ter models to cover the range in possible initial conditions. 
We present the results from a number of these simulations. 



listed in Table |^, in order to illustrate the relevant physical 
processes of accretion in stellar clusters. The simulations are 
characterised by the number of stars, the gas fraction, the 
stellar distribution and the number of Jeans masses in the 
gas, 

Mj=[^-^\' {\^P)'"\ (1) 



V 2Gn 



where Rg is the gas constant, G is the gravitational con- 
stant, and ^ is the mean molecular weight. Other cluster 
parameters are the initial mean Mach number of the stars 
and the distribution of stellar masses. 

All the clusters considered here are initially spherical 
and have uniform gas distribution. In all cases, the gas is 
initially static and has only thermal support. The stars ini- 
tially hav a small virial ratio in all cases. The initial condi- 
tions are termed 'cold' when they contain at least as many 
Jeans masses as stars. Such initial conditions are the most 
likely if the cluster forms through a fragmentation event. An 
isothermal gas equation of state is used throughout. This 
means that, in the absence of accretion, the gas will collapse 
to a singularity in approximately one free-fall time. 

The initial spatial distribution of stars is either uniform 
or n{t) oc as might be expected after a violent relaxation. 
One simulation (Run C) was performed with an initial range 
of stellar masses in order to ascertain the importance of the 
initial stellar mass distribution. 



3 ACCRETION AND CLUSTER DYNAMICS 

One of the first things we need to understand in terms of the 
dynamics of accretion in clusters is the timescale. This sets 
the mean accretion rate and thus how important accretion 
can be in setting the final stellar mass. Gas accretion by a 
star is given by the general formula 



(2) 



where p is the gas density and v^^i is the relative gas-star 
velocity and iiacc is the accretion radius. 

One possible estimate of the accretion timescale is using 
Bondi-Hoyle accretion (Bondi & Hoyle 1944; Bondi 1952), 
corresponding to an isolated star accreting from a uniform, 
non-self-gravitating medium. In this model, the accretion 
radius is given by 

CbhGM, 



+ ci 



(3) 



where M, is the stellar mass, and Cs is the gas sound speed 
and the constant Cbh is generally taken to be Cbh ~ 2 
(Bondi & Hoyle 1944). This approach neglects the self- 
gravity of the gas, the presence of other stars, the cluster 
potential and how these affect the accretion. Estimating the 
accretion timescale as 



the Bondi-Hoyle formalism implies 
If the cluster is virialised 



(4) 



(5) 
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Run 


No. of stars 


Gas fraction 


stellar dist" 


N(M,j) 


Mach no. 


Initial masses 


Run A 


100 


90 


uniform 


120 


0.5 


equal 


Run B 


100 


90 




2.8 


0.2 


equal 


Run C 


100 


90 


uniform 


120 


0.5 


range 


Run D 


1 


99 




130 


0.2 




Run E 


1 


99 




130 


2. 




Run F 


30 


82 


uniform 


110 


0.6 


equal 


Run G 


100 


90 


r-2 


120 


0.5 


equal 



Table 1. The initial cluster models used to illustrate the physics of the accretion process are listed by the number of stars in the cluster, 
the percentage of total mass initially in gas, the stellar distribution, the number of Jeans masses contained in gas, the mean Mach number 
of the stars and the distribution of initial stellar masses. 



-Rclust 

then 



tacc{BH) ^ Nt. 



(6) 



(7) 



where N is the number of stars in a cluster and tdyn is the 
cluster's dynamical or crossing time. In studies of accretion 
in small-N clusters (Bonnell et al. 1997), this timescale was 
found to be much too long as it neglects the motion of the 
gas in the cluster. 

Another possible estimate is obtained by considering 
how the gas evolves under the cluster potential. Assuming 
the gas is gravitationally unstable (otherwise it would not 
have formed stars in the first place), dominates the cluster 
potential, lacks rotational support and remains isothermal, 
the gas falls towards the centre of the cluster on the free-fall 
timescale. 



iff = 



32Gp' 



(8) 



In the absence of any accretion onto existing stars, the gas 
collapses and forms an additional star containing all the gas 
mass in approximately one free-fall time. This sets an up- 
per limit on the accretion timescale as the gas remains dis- 
tributed in the cluster for no more than approximately one 
free-fall time. 

During the collapse, the gas density increases and the 
cluster potential becomes deeper. Thus as long as the stars 
remain inside the gas sphere, the increase in the gas density 
will increase the accretion rate (equation 2) . As the density 
increases towards infinity on a dynamical timescale, the ac- 
cretion rate also increases on this timescale till all the gas 
is accreted (or has formed another star at the cluster cen- 
tre). Thus, the accretion timescale in a self-gravitating gas 
is given by 



ts- 



(9) 



Comparing these estimates to the numerical simula- 
tions, we find that the accretion occurs on the dynamical 
timescale as the gas density increases dramatically on this 
timescale. Figure ^ plots the evolution of the total accreted 
mass versus time for Run A (a cold cluster containing 100 
stars with an initial gas fraction of 90 per cent). The stars 
initially have a total mass of 0.1. They accrete somewhat 
slowly over the first half of a free-fall time and only double 
their mass by 0.8 tg. By 0.95 tfi the stars have accreted all 
the gas in the cluster and have increased their mean mass 
tenfold. 
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Figure 1. The total accreted mass (solid line) is plotted as 
a function of time for Run A (a 'cold' cluster containing 100 
stars and 90 per cent of its initial mass in the form of gas) . The 
half-mass radius of the stars is also plotted (dashed line). Note 
that both the gas accretion and the cluster evolution occur on a 
dynamical timescale. 



The half-mass radius of the stellar distribution is also 
plotted in Figure |l| showing how the cluster evolves on a 
dynamical timescale. The stars are initially cold and thus 
in the absence of any gas accretion would revirialise at half 
their initial radius. Instead, the stars' dynamical evolution 
is due to the combination of the accretion and the changing 
cluster potential due to the collapse of the gas. As the gas 
collapses, the deepening of the cluster potential forces the 
stars to follow the gas and shrink. Additionally, the accre- 
tion increases the stellar masses making them more bound 
to each other which causes the stellar distribution to shrink 
further. Thus not only does the accretion occur on a dy- 
namical timescale but the cluster as a whole evolves on this 
timescale. This evolution may be important in forming mas- 
sive stars if they cannot accrete, due to their large radiation 
pressure, past « IOMq. In this case, the cluster shrinkage 
due to the accretion can be sufficient to induce stellar col- 
lisions which form the most massive stars (Bonnell, Bate & 
Zinnecker 1998). Stellar collisions are not included in the 
simulationss presented here as the finite stellar size is not 
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Figure 2. The accretion rate onto 100 stars in a cluster is plotted 
as a function of the radius in the cluster of each of the stars 
(Run B). The accretion rates are averaged over O.ltfj and are 
calculated at t fx 0.65%. The accretion rate increases towards 
the centre of the cluster which leads to mass segregation. Initially 
the cluster is centrally condensed and contains 10 per cent of its 
mass in stars. The accretion rate is in units of the total cluster 
mass per tg. 



Figure 3. The resultant stellar masses (in units of the total 
cluster mass) after accretion are plotted as a function of each 
star's position in the cluster at t fa Itg for Run A. The initial 
stellar mass is = 0.001. The mass segregation is apparent 
even though there are significant numbers of low-mass stars near 
the cluster centre. The cluster is initially uniform and cold and 
contains 10 per cent of its mass in stars. 



modelled. The reader is referred to Bonnell et al. (1998) 
for a discussion on how accretion-induced shrinkage of the 
clsuter can result in stellar collisions. 

Although many of the clusters investigated here are cold 
and contain ~ 100 Jeans masses, they do not in general col- 
lapse to form additional stars. Instead, the gas was accreted 
onto individual stars that the gas encountered on its way 
to the cluster centre. Exceptions to this were the simula- 
tions run with only one star present in the gas cloud. In 
these cases the gas can easily collapse to the centre without 
encountering (being accreted by) the single star. 



4 ACCRETION AND MASS SEGREGATION 

One of the goals of studying accretion in stellar clusters is 
to determine if it can explain the initial mass segregation 
that is observed in young stellar clusters (e.g. Hillenbrand & 
Hartmann 1998). This segregation of the more massive stars 
towards the centre of the cluster is not due to two-body re- 
laxation in the cluster as the clusters are often too young 
(Bonnell & Davies 1998). Moreover, simple Jeans mass ar- 
guments imply that the lowest mass stars should be found 
near the centre of the cluster where the gas density is high- 
est, not the most-massive ones (Zinnecker, McCaughrean & 
Wilking 1993; Bonnell et al. 1998). Simulations of accretion 
in small stellar clusters (Bonnell et al. 1997) found that the 
accretion rate was generally higher for the stars nearer the 
cluster centre. In this section, we investigate how the accre- 
tion rate depends on the stellar positions in the cluster and 
how this can result in a mass-segregated cluster. 

Accretion leads to mass segregation in two ways. Firstly, 



as the system relaxes to a centrally-condensed distribution, 
the gas density increases towards the centre of the cluster 
which then results in higher accretion rates in the centre 
than elsewhere. Secondly, mass accretion onto an individual 
star decreases its kinetic energy (mass loading) which forces 
it to sink deeper into the cluster potential. Thus, stars that 
accrete gas sink towards the centre and stars that are in the 
centre accrete more gas as the density is highest there. 

Figure ^ plots the individual accretion rates versus stel- 
lar position in the cluster for Run B. These accretion rates 
are averaged over approximately O.ltg at time t ~ 0.65tg. 
The 100 stars are initially centrally condensed in this model 
and the gas, comprising 90 per cent of the total mass, is rel- 
atively warm in that it contains only 2.8 Jeans masses. The 
stars dominate the central part of the potential due to their 
condensed distribution. Figure ^ shows that the accretion 
rate increases with decreasing radius with those stars near 
the cluster centre accreting at a rate more than a hundred 
times those near the outside. This differential accretion is 
due to the higher gas densities found near the centre of the 
cluster and obviously leads to a mass-segregated cluster on 
the accretion (dynamical) timescale. The most massive stars 
are the ones that started out nearest the cluster centre. They 
have the largest accretion rates and will therefore always be 
nearest the cluster centre (unless ejected by a dynamical 
interaction) . 

For clusters whose stars are initially uniformly spa- 
tioally distributed, the gas density is initially uniform as well 
and only becomes significantly centrally concentrated on a 
free-fall timescale. In this situation, it is not necessarily the 
stars initially nearest the cluster centre which accrete the 
most. Instead, it is those which, regardless of how close they 



© 0000 RAS, MNRAS 000, 000-000 



Accretion in stellar clusters 5 



0.1 



m 
m 



0.01 



10-3 



10- 




10-3 



mass 



initial 



Figure 4. The intermediate (i fa 0.87tff, open pentagons) and 
final (t !^ 1-Otfi, filled triangles) masses are plotted versus the 
initial stellar mass for Run C, a cluster containing 100 stars with 
a distribution of initial masses. The solid line denotes the initial 
range of masses. The cluster is 'cold' and contains 90 per cent of 
its initial mass in the form of gas. 



are to the cluster centre, are in a favourable location (e.g. 
less local competition) that accrete the most. These stars 
then sink towards the centre and continue to accrete at a 
higher rate. This still results in a mass segregated cluster 
but there is no direct correlation between the stars' initial 
positions in the cluster and their eventual mass. An example 
of the resulting mass distribution in the cluster is given in 
Figure |^ for Run A (an initially uniform cluster containing 
100 stars with 90 per cent of its mass in the form of cold 
gas and containing 120 Jeans masses). We note here that al- 
though the trend of higher mass stars towards the centre is 
very distinct, there are still low- mass stars near the cluster 
centre as is observed in young stellar clusters such as the 
ONC (Hillenbrand 1997; Hillenbrand & Carpenter 2000). 



5 ACCRETION AND STELLAR MASSES 

We have shown above that accretion can result in a spectrum 
of stellar masses and in mass segregation even when the 
initial mass distribution is that of equal masses. A further 
question is how dependent the final masses are on the initial 
mass distribution. We performed a test case with a spread 
of a factor of five in the initial masses. Figure ^ shows the 
relationship between the initial and post-accretion stellar 
masses for Run C, an initially uniform, 'cold' cluster with 
90 per cent of its mass initially in the form of gas. The post- 
accretion mass distribution is shown when the stars have 
increased their total mass (from 10 per cent) to 25 percent 
and 100 per cent of the cluster mass. We see that the initial 
masses play a large role in determining the masses at the 
intermediate stage, ie that there is still a strong correlation 
between initial and post accreted masses. On the contrary, 
at the later stage where the stars have accreted all the gas. 



the correlation between final and initial masses is much less 
clear. There is still a general trend that the mean final mass 
increases with the initial mass, but the scatter at any given 
initial mass is so large (a factor 20) that a star's final mass 
places only the weakest constraints on the initial mass of the 
condensation from which it formed. 

The reason for this difference in the intermediate and 
final mass distribution lies in the cluster evolution. The clus- 
ter is initially uniform (constant density) such that the only 
difference in the accretion rates comes from the accretion ra- 
dius. This is dependent on the stellar masses (see below) and 
as such leads to a dependence of the post-accretion mass on 
the initial mass. The cluster is unstable in a uniform configu- 
ration and evolves (in this case due to gravitational collapse) 
to a more centrally-condensed configuration. Once the clus- 
ter is centrally condensed, the gas density play a more im- 
portant role in the accretion rate such that the mass accreted 
is more due to the gas density than due to the accretion 
radius. This is especially so if the accretion radius is only 
weakly dependent on the star's mass as we shall see below. 
Thus, a centrally-concentrated cluster allows even initially 
lower-mass stars to have higher accretion rates, and eventu- 
ally attain higher masses if they are in regions with high gas 
density. 



6 AN ANALYTIC PRESCRIPTION FOR 
COMPETITIVE ACCRETION 

In addition to seeing how accretion results in mass segre- 
gation and a mass spectrum (Figure ^), one of our goals 
is to understand the physics of the accretion process. This 
involves determining what sets the accretion rates of indi- 
vidual stars. A parametrised form for the accretion rates will 
allow the accretion process to be applied to larger groups of 
stars without needing to resolve the accretion fiow around 
each star. 

The first possible parametrisation of the accretion rate 
is the Bondi-Hoyle formalism discussed above. In this case, 
the accretion radius is basically the radius where the gravita- 
tional energy due to the star is larger than the kinetic energy 
( including the unperturbed gas velocity relative to the star, 
iirci, and sound speed, Cs). Thus (Bondi & Hoyle 1944), 



2GM, 



and 



(10) 



(11) 



Alternatively, the cluster potential sets another radius, the 
tidal radius, due to the multiple gravitational sources. This 
radius is analogous to the Roche-lobe radius in accreting 
binary systems. We adopt this radius as our second potential 
accretion radius. This tidal-lobe radius is 



-n-tidal — (-^tidal i -77 I ' 



(12) 



where M« is the star's mass, and A/cnc is the mass enclosed 
within the cluster at the star's position r«. In the Roche- 
lobe approximation, Ctidai ~ 0.5 (e.g. Pacynski 1971) which 
we will adopt throughout this paper. The accretion rate is 
then 
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Figure 5. The stellar velocity dispersion (open pentagons) and 
the mean relative gas velocity in the tidal-lobes of the individual 
stars (filled triangles) are plotted in units of the gas sound speed 
Cs as a function of time for Run A . The cluster contains 100 
stars, is initially uniform (stars and gas) and contains 90 per cent 
of its mass in gas (120 Jeans masses). 

M* « 1^ pVralRlidal- (13) 

The expectation is that the smaller of the two radii 
should determine the accretion rate as it will be this radius 
which ultimately decides if the gas is bound to the star and 
thus will be accreted. Considering the high stellar velocities 
in a cluster, the Bondi-Hoyle radius is commonly expected 
to be the smaller of the two. 

One difficulty is estimating the relative velocity which 
goes into both the Bondi-Hoyle radius (equation 10) and 
the Bondi-Hoyle accretion rate (equation 11). We make two 
estimates of this velocity, one calculated as the velocity of 
each star relative to the centre of mass of the cluster and 
the second is the velocity relative to the gas at the tidal- 
lobe radius. This relative velocity is calculated by averaging 
each component of the three-dimensional velocity vectors of 
the particles which overlap the tidal radius. We shall refer 
to the Bondi-Hoyle accretion estimate using this relative ve- 
locity as the modified Bondi-Hoyle accretion. The evolution 
of each of these velocities is shown in Figure ^ for Run A, a 
cold cluster of 100 stars initially comprising 10 per cent of 
the total mass. Both velocities are initially subsonic but the 
stellar velocities quickly become supersonic as the cluster 
collapses. In contrast, the relative gas velocity at the tidal- 
lobe remains subsonic throughout the evolution. This occurs 
as the gas is accelerated under the same potential that the 
stars are. Thus, the gas tracks the stellar motions and the 
relative velocity is small. 

The gas density is also estimated at the tidal-lobe radius 
as that is the radius at which the density is unaffected by 
the star's potential. This density is calculated as the average 
over the particles which overlap the tidal radius. Estimates 
of the gas density and relative velocity at the tidal-radius are 
generally well determined as the number of SPH particles 




0.2 0.4 0.6 0.8 1 
time (tjj) 

Figure 6. The mean accretion rate (filled circles) is plotted 
as a function of time for Run A, a cluster of 100 stars which 
initially has 90 per cent of its mass in gas and is cold. Estimates 
of the accretion using tidal-lobe as the accretion radius (filled 
triangles) and using the Bondi-Hoyle radius (open pentagons) are 
also plotted. The upper Bondi-Hoyle curve uses the relative gas 
velocity in each star's tidal lobe (modified Bondi-Hoyle) whereas 
the bottom curve uses each star's velocity in the cluster. 



there is typically large. Exceptions to this occur when a 
star is ejected from the cluster and the gas density falls 
to zero. Another problem that was encountered was due to 
the presence of other stars in the tidal-lobes. The tidal-lobe 
formulation does not allow for the effects of close binary 
systems (with i?scp < Rtnn^i) so that the gas velocities and 
density do not correspond to only one star. Such systems are 
generally excluded from the analysis that follows (except in 
§ 6.3). 

Using the above estimates of the accretion rate as a 
function of the local gas density, relative stellar velocity and 
the accretion radius, we can compare these with the accre- 
tion rate determined from the SPH simulations. Figure ^ 
plots the evolution of the mean accretion rate for Run A, 
a 'cold' cluster of 100 stars where 90 per cent of the initial 
cluster mass is in the form of gas. Only those stars not in bi- 
naries are included. The comparison includes two estimates 
of the Bondi-Hoyle radius using each star's velocity relative 
to the cluster's centre of mass and its velocity relative to the 
gas in the tidal-lobe. 

The mean accretion rate determined by the SPH simu- 
lations increases with time until near the end of the simula- 
tion when the gas is significantly depleted. The Bondi-Hoyle 
accretion rates are both initially much higher than the SPH 
accretion rate as the stars are initially at rest and the gas 
is cold. The modified Bondi-Hoyle accretion rate (using the 
velocity relative to the local gas) remains much higher than 
the SPH accretion rate. This is because, as shown in Fig- 
ure ^, this velocity remains subsonic throughout the simula- 
tion, and thus Rbh > RtiAai- On the other hand, the Bondi- 
Hoyle accretion rate using the stellar velocity relative to the 



© 0000 RAS, MNRAS 000, 000-000 



Accretion in stellar clusters 7 



100 




0.1 - 



0.01 E- 



0.001 E- 



radius 



10 E 



1 E 



T3 



0.1 



10- 



10 



-3 _ 




-0.5 0.5 
projected distance 



Figure 7. The accretion rate (filled circles) of a single star in 
a cloud of cold gas is plotted as a function of its radius in the 
cloud (Run D). The star initially contains one percent of the 
total mass. Estimates of the accretion rate based on the tidal- 
lobe (filled triangles) and modified Bondi-Hoyle (stars) are also 
plotted. The open squares denote the tidal-lobe accretion that 
includes a component of spherically-symmetric infall. 

cluster centre of mass decreases as the stars accelerate and 
move supersonically. The accretion rate is thereafter signifi- 
cantly smaller than the SPH accretion rate. In contrast, the 
tidal-lobe accretion rate is generally very similar to the SPH 
accretion rate throughout the simulation. 

From this comparison we can see that the Bondi-Hoyle 
accretion rate is either too high or too low depending on 
which velocity is used. It is also readily apparent why the 
naive accretion rate discussed in Section 3 gives a timescale 
that is too long. The stellar velocity relative to the clus- 
ter's centre of mass is not the appropriate velocity. It is the 
velocity relative to the local gas from which the star is ac- 
creting which must be used. Thus, from now on, we shall 
discuss only the modified Bondi-Hoyle accretion rate using 
the stellar velocity relative to the local gas. 

6.1 Single-star clusters 

In order to ensure that we can correctly model the accre- 
tion rate, we need to compare the different parametrisations 
of the accretion for individual stars. Figure ^ plots the ac- 
cretion rate onto one star in a cold gas cloud where 99 per 
cent of the initial mass is in the form of gas (Run D). Thus, 
the gas provides the potential as would do an entire cluster 
but without the difficulties of binary systems and also with 
increased numerical resolution. The accretion rate onto the 
star is plotted as a function of it's position (radius) in the 
gas cloud. The accretion rate increases as the star falls in 
towards the centre of the cloud, due to the increased density 
there as the cloud collapses to a centrally-condensed config- 
uration. Also plotted in Figure ^ are the modified Bondi- 
Hoyle accretion rate (using the velocity relative to the gas) 



Figure 8. The density of SPH gas particles around a single-star 
cluster (Run E) are plotted as a function of the projected distance 
relative to the direction of the star's motion. The star is travelling 
inwards at Mach 0.9 relative to the local gas. The size of the tidal 
lobe is indicated and corresponds approximately to where the 
density starts to increase towards the star. The modified Bondi- 
Hoyle radius is larger than the plotted region (/Jbh ~ 2.). Note 
that the gas density is higher at projected distances < 0, indicat- 
ing that the gas is predominantly accreted from behind the star. 
The slight density gradient is because the star is infalling a bit 
faster than the gas and thus seeing the higher density ahead. 

and two estimates of the tidal-lobe accretion rate. The first 
of these estimates includes only the relative gas velocity, 
Vy-ci, whereas the second (and higher) estimate includes an 
adhoc estimate of the degree of spherical infall through the 
tidal- lobe, t;inf. This spherical accretion can occur over 47r 
sterradians, but, as a first approximation, we use and in- 
termediate value of Stt as the accretion is not completely 
spherical. Thus 

M, = 7rpiiLc('^rcl + 3Uinf), (14) 

where Winf is the mean infall velocity of particles at the tidal- 
lobe radius, Winf = v • r/|r|. 

As seen above, the modified Bondi-Hoyle accretion rate 
is much larger than the SPH accretion rate. The tidal-lobe 
estimates are much closer to the numerically determined 
value. The estimate not including spherical infall is slightly 
lower than the SPH value whereas the estimate including 
the spherical accretion is very slightly higher than the SPH 
value. 

The tidal-lobe accretion is a good fit to the accretion 
rate even in simulations where the star is initially super- 
sonic relative to the cluster's centre of mass. This occurs 
due to the convergence of the stellar and gas velocities as 
the cluster collapses. Both the star's and the gas velocities 
are increasingly dominated by the gravitational acceleration 
of the cluster. As both the star and the gas start collaps- 
ing towards the centre of the system, the gas velocity in- 
side the star's tidal-lobe becomes more similar to the star's 
velocity. Thus, the relative velocity decreases and the fiow 
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Figure 9. Estimates of the accretion rate using the modified 
Bondi-Hoyle (open pentagons) and tidal-lobe (filled triangles) for- 
malisms are plotted against the SPH determined accretion rates 
for each star in a cluster of 30 (Run F). Only stars that do not 
have companions within -Rtidal S'i's included. The cluster is ini- 
tially uniform, cold and contains 82 per cent of its mass in gas. 

becomes subsonic. This effect is aided by the accretion as it 
helps converge the star's velocity to that of the flow. Thus, 
RtiAsi < Rbh and the tidal-lobe accretion determines the 
accretion rate. 

The finding that the tidal-lobe radius is a good approx- 
imation to the accretion radius is further backed up by look- 
ing at the gas density near the accreting star. Figure |8| plots 
the gas density as a function of the projected distance from 
the star for a single-star initially travelling at Mach 2 in 
a cluster potential dominated by gas (Run E). Figure ^ is 
plotted at a time by which the gravitational acceleration of 
the cluster has converged the stellar and gas motions such 
that both the star and the gas are moving inwards and the 
relative velocity of the local gas is Mach 0.9. The projected 
distance is calculated as the distance relative to the star's 
velocity, 

dproj = (15) 

We see that the gas density starts to increase at a distance 
that corresponds roughly to the tidal-lobe radius as expected 
if this is the physical accretion radius. It is also worth noting 
that the gas density is somewhat skewed relative to the star's 
position with higher density material at greater distances 
behind the star then in front of it (relative to its motion). 
This occurs as generally the material is accreted from behind 
the star, as in the classical Bondi-Hoyle accretion. 

6.2 Many-star clusters 

Accretion onto stellar clusters is considerably more complex 
than the above models of a single star in a gas-dominated 
potential. In a stellar cluster, accretion is not only limited 



by the tidal potential of the overall cluster but also by com- 
petitive accretion from other stars. Figure ^ plots the com- 
parison of the SPH accretion rate with that based upon 
tidal-lobe and modified Bondi-Hoyle accretion for Run F, a 
'cold' cluster of 30 stars initially containing 82 per cent of 
its mass in gas. The accretion rates are plotted for all stars 
that do not have another star within their tidal-lobe. This is 
done as the tidal-lobe determination neglects the presence 
of any companions (the effect of other stars in the tidal-lobe 
is discussed in §6.3). 

The evolution of the accretion rate for the cluster is 
from low accretion rates towards higher accretion rates as 
the cluster collapses under its self-gravity. The quantisation 
of the lowest accretion rates is due to numerical accuracy 
and therefore no conclusions are based upon these points. 
We see that, as in the case of the single-star cluster, the 
tidal-lobe accretion is a good approximation to the numeri- 
cal accretion rate determined by the SPH code. In contrast, 
the modified Bondi-Hoyle accretion rate, using the relative 
velocity inside the tidal-lobe, is much too high. It is impor- 
tant to note that the tidal-lobe accretion is a good approx- 
imation to the SPH determined accretion rate, although it 
generally slightly underestimates it. In contrast, not only is 
the modified Bondi-Hoyle accretion rate too high, it also di- 
verges from the numerical accretion rate towards the higher 
accretion rates. The tidal-lobe accretion is a better fit to the 
accretion rate as -Rtidai < -Rbh when the local gas velocity 
is used to calculate Rbh- 

This agreement between the numerical and tidal-lobe 
accretion rates was found for most clusters as long as the 
gas is reasonably cold. If the gas is sufficiently warm then 
Rbh < -Rtidai and the Bondi-Hoyle is a reasonable fit to the 
accretion. Such a situation occurs if the gas contains only 
a few Jeans masses and thus should not have been able to 
fragment to form the number of stars found in the cluster. 
One possibility for a cluster with relatively warm gas is when 
feedback from young stars heats up the gas. Other excep- 
tions were found when many stars are in the same tidal-lobe 
such as occurs when gas infalls from outside the cluster onto 
a star-dominated potential, and virialised cluster (see be- 
low). Generally up to a third of systems have an other star 
within their tidal-lobes at any given time, although this can 
be significantly higher in the stellar dominated potentials. In 
this last case, most of the stars have -Rbh < J2tidai although 
some are still better modelled by the tidal-lobe accretion 
when their velocities are inwards. 

6.3 Accretion in a stellar dominated potential 

As noted above, companions within a star's tidal-lobe com- 
plicate the competitive accretion picture. This becomes in- 
creasingly important during the evolution as the gas is ac- 
creted and the stars come, at least in the centre, to domi- 
nate the cluster potential. In such circumstances, compan- 
ion stars in each other's tidal-lobes becomes common. In 
such circumstances the tidal-lobe accretion formalism as de- 
scribed above does not work as the mass crossing into an 
individual tidal-lobe can still be accreted by several objects. 

When a star has one or more companions inside its nom- 
inal tidal-lobe, it's velocity can be significantly higher than 
the surrounding gas as it is under the additional acceleration 
of the relatively close companion(s). The extreme of this is 
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Figure 10. The ratio of the theoretical accretion rates compared 
to that calculated by the SPH code is plotted against the ratio 
of the tidal-lobe radius to the modified Bondi-Hoyle radius for 
Run C. The top panel uses the modified Bondi-Hoyle accretion 
radius whereas the bottom panel uses the tidal-lobe radius as the 
accretion radius. The middle panel uses the tidal-lobe accretion 
radius but adapted to include the eff'ects of multiple stars inside 
this radius (see text). The tidal-lobe radius gives a better fit to 
the accretion rate when Rtidal 5; ^BH- The accretion rates are 
calculated for a cluster of 100 stars with variable initial masses 
and an initial gas fraction of 90 per cent. 

when the gas is falling onto a virialised cluster and hence the 
stellar and gas velocities are generally uncorrelated. In this 
case, the modified Bondi-Hoyle radius becomes smaller than 
the tidal-lobe radius {Rbh < Rudai) and thus determines the 
accretion rate. Figure |l^ plots the ratio of the theoretical to 
numerical accretion rates versus the ratio of the tidal-lobe 
radius to the modified Bondi-Hoyle radius for Run C. The 
tidal-lobe accretion (bottom panel) models well the accre- 
tion rate when the tidal-lobe is smaller than the modified 
Bondi-Hoyle radius (iitidai/^BH ^ 1). When the tidal-lobe is 
significantly larger than the modified Bondi-Hoyle radius, as 
is the case when other stars are present inside the tidal-lobe, 
then the tidal-lobe formalism seriously overestimates the ac- 
cretion rate. The top panel of Figure hd shows the case when 
modified Bondi-Hoyle accretion is used. As found above, the 
modified Bondi-Hoyle formalism seriously overestimates the 
accretion rate when i?tidai/J?BH i$ 1 as the accretion is deter- 
mined by the size of the tidal-lobe. In contrast, the modified 
Bondi-Hoyle formalism works better than the tidal-lobe for- 
malism when -Rtidai >> Rbh- The scatter towards the right 
of each panel in Figure ^ is due to numerical resolution 
(low accretion rates) and sampling as the effective accretion 
radius becomes small (few particles to calculate Urci , p) ■ 

The middle panel of Figure ^ shows the result when 
tidal-lobe accretion is adapted to include the effects of com- 
panions. This is done by estimating the fraction of the ma- 
terial enterring the tidal-lobe which is accreted by the star. 
In the tidal-lobe formalism, this can be calculated by mod- 
ifying equation (13) as 



Figure 11. The Mach number of the relative gas flow near each 
star is plotted as a function of the stars' positions in the cluster for 
Run G (lower panel). The upper panel shows the corresponding 
ratio of the tidal to modified Bondi-Hoyle radius as a function of 
position in the cluster. The stars dominate the inner parts of the 
cluster such that the stellar and gas velocities are uncorrelated, 
resulting in iJtidal > ^BH- The gas dominates the outer parts of 
the cluster such that the gas and stellar velocities are correlated 
and -Rtidai < Rbh- The values are averaged over O.ltff at t 
l.Otff. The relative overdensity of stars a.t r 0.12 indicates that 
the deepest part of the cluster potential, where the most massive 
stars are, is not at the centre of mass of the system. 



/ MJ' \ 2 

M, = TVpVini 2- -Rtidai (16) 

where Vh-ii is the mean spherical infall through the tidal-lobe, 
and A'^comp is the number of stars in -Rtidai- The mass term 
effectively recalculates the cross-section of the tidal-lobe to 
include the companions (as -Rtidai oc (M,/Menc)'''^), while 
using Vint removes the effect of the orbital velocity as the 
stars move around each other inside the tidal-lobe. We see 
this estimate of the tidal-lobe accretion works fairly well for 
J^tidai S; 10-Rbh, and models to first order the reduction in 
the accretion rates due to companions inside the tidal-lobe. 
The limitation of this approach is that it neglects the indi- 
vidual velocities of the stars which will modify the accretion 
rates. Thus, once -Rtidai >> Rbh we are left with using the 
modified Bondi-Hoyle formalism. 

As a stellar cluster accretes from the gas, the mass- 
fraction in the stars increases while it decreases for the gas. 
Once the stars dominate the cluster potential, their evolu- 
tion is no longer determined primarily by the collapsing gas. 
They then virialise and thus have velocities uncorrelated to 
the infalling gas. This transition is illustrated in Figure |ll| 
showing the relative gas velocity as a function of stellar po- 
sition for Run G. The stars are initially centrally condensed 
which results in their dominating the central portions of 
the cluster potential. In the outer regions of the cluster, the 
'cold' gas still dominates the potential such that both stars 
and gas are infalling and the relative gas velocity is subsonic. 
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Inside of r = 0.9, the gas contributes less than 25 per cent 
of the mass and thus no longer determines the dynamics. 
The stars become virialised and the relative gas velocities 
are large as the stellar and gas motions are uncorrelated. 
Figure nil also plots the ratio of the tidal-lobe radius to the 
modified Bondi-Hoyle radius as a function of position in the 
cluster. Stars exterior to r = 0.9 have Rtidai < Rbk as the 
relative gas velocities are subsonic whereas stars interior to 
r = 0.9 have -Rtidai > Rbh as the relative gas velocities are 
supersonic. Exceptions to this are the most massive stars in 
the centre of the cluster which, due to their high mass, have 
low velocities and thus -Rtidai ^ -Rbh once again. 

In general as a cluster accretes from its gas reservoir, 
the stars begin to dominate the potential from the inside 
out. This occurs due to the higher accretion rates in the 
centre and due to the settling of the accreting stars. Thus, 
as the gas infalls from larger radii onto the cluster it will pass 
through two regimes. First, where the gas dominates the po- 
tential, and the stellar dynamics, the accretion is dominated 
by the tidal-lobes of each star. Second, further in where the 
stars dominate the potential and are virialised, the accretion 
is determined by a Bondi-Hoyle accretion. 



7 DISCUSSION OF ACCRETION MODELLING 

In the above sections, we found that gas accretion in stellar 
clusters is generally well modelled by using the tidal-lobe 
as the accretion radius. This formalism works as the clus- 
ter dynamics are dominated by the overall collapse of the 
cluster due to its self-gravity. In this case the gas in the 
tidal lobe of a star has a similar velocity to the star as 
both are primarily due to the same gravitational acceler- 
ation. This results in a large modified Bondi-Hoyle radius 
such that the smaller tidal-lobe radius dominates the accre- 
tion, 7?acc ~ iitidai < Rbh- Such a sceuario will always result 
when the gas dominates the gravitational potential and is 
dynamically moving under the influence of its self-gravity. 
In such a case, even if the stars are initially virialised, the 
changing potential and the gas accretion will maintain sim- 
ilar gas and stellar velocities. 

Complications arise when more than one star is in a 
tidal-lobe region as this increases the stellar velocity rela- 
tive to the gas velocity in the tidal-lobe. Gas accretion can 
then be modelled by including the additional stars in the 
calculation of the tidal-lobe (and thus the accretion cross 
section). Once the stars dominate the gravitational poten- 
tial and are virialised, then -Rtidai >> ^^bh and accretion of 
infalling gas can be well modelled with modified Bondi-Hoyle 
accretion although some care has to be taken to calculate 
the appropriate density and velocity of the gas. 



8 CONCLUSIONS 

We find that accretion in stellar clusters occurs on a dynami- 
cal timescale. Both the gas and stellar motions are governed 
by the cluster's gravitational potential. This ensures that 
the gas will fill the depleted areas and thus find the stars 
onto which it accretes. 

In clusters where gas dominates the potential, the ac- 
cretion process is better modelled by using the tidal-lobe 



radius, -Rtidai, as the accretion radius rather than the com- 
monly used Bondi-Hoyle radius, Rbh- This works as both 
the stellar and gas velocities are dominated by the the gas- 
dominated gravitational potential, and thus have similar ve- 
locities. In this case, we have Titidai ^ -Rbh and the smaller 
^tidai determines the accretion rate. This will occur as long 
as the gas is the dominant component and is free to col- 
lapse (is unsupported against gravity). In contrast, where 
the stars dominate the cluster potential and are virialised, 
the stellar and gas velocities are uncorrelated resulting in 
high relative gas velocities and -Rbh < Rtidai- In this case 
the accretion is better modelled by Bondi-Hoyle accretion. 

Accretion in a stellar cluster is highly non-uniform with 
different stars accreting at significantly different rates. Cou- 
pling this with a large mass fraction in the form of gas in 
the cluster results in a spectrum of stellar masses, even from 
initially equal masses. The accretion rate depends primar- 
ily on the local gas density and thus on the star's position 
in the cluster as the gas quickly concentrates towards the 
deepest part of the potential well. There is a weak depen- 
dence on the mass (.M* cx M'^^^ for tidal-lobe accretion) such 
that lower initial masses bias against accretion. This can be 
overwhelmed by higher gas densities found in the centre of 
the cluster. Thus, competitive accretion naturally results in 
a mass-segregated cluster on the formation timescale and 
does not require subsequent two-body relaxation. 
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